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A  basic  program  is  presented  to  allow’  high  speed  digital  compu¬ 
tation  of  continuous  transformation  diagrams  from  data  available  in 
the  form  of  isothermal  transformation  diagrams.  The  program  follows 
the  method  of  Grange  and  Kiefer  and  continuous  transformation  diagrams 
computed  by  this  method  are  shown  to  agree  favorably  with  experimental 
diagrams. 
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Much  of  the  available  information  on  which  the  heat  treatments  of 
a  particular  steel  are  based  is  found  in  the  isothermal  transformation 
diagram.  Such  is  the  case  despite  the  fact  that  the  usual  thermal  treat¬ 
ments  involve  continuous  transformations,  i.e.  transformations  which  occur 
on  cooling,  rather  than  isothermal  transformations.  This  rather  unusual 

practice  is  continued  because  of  the  abundance  of  published  isothermal 

(1  2) 

transformation  diagrams  9  (which  are  relatively  easy  to  determine  ex¬ 
perimentally)  and  the  difficulty  of  experimentally  determining  the  con- 

(3  4) 

tinuous  transformation  curves  9  ,  but  can  easily  lead  to  the  production 

of  improper  microstructures  and  may  not  allow  full  advantage  to  be  taken 
of  the  inherent  hardenability  of  a  steel. 

A  mathematical  method  of  estimating  the  continuous  transformation 
diagram  from  the  available  isothermal  data  was  developed  by  Grange  and 

(5) 

Kiefer  ;  but,  due  to  the  mathematical  approximations  required  in  the 
calculations,  the  computation  of  continuous  transformation  diagrams  has 
been  limited  to  individual  cases  and  no  large  scale  compilation  of  these 
diagrams  is  available. 

With  the  present  availability  of  high  speed  digital  computers,  the 
problem  of  time  consuming  approximations  may  be  minimized  and  continuous 
transformation  diagrams  may  be  constructed  with  the  insertion  of  isothermal 
data  into  a  basic  program.  It  is  the  purpose  of  this  report  to  announce 
the  availability  of  such  a  program  and  compare  the  computed  and  experimental 
continuous  transformation  diagrams  for  two  different  steels. 

(5) 

Let  us  first  review  the  method  of  Grange  and  Kiefer  and  then  compare 
the  results  of  the  machine  calculations  with  the  experimentally  determined 
diagrams.  In  order  to  apply  this  method,  one  requires  the  isothermal 
transformation  diagram  and  a  set  of  curves  representing  cooling  at  constant 
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rates.  A  section  of  an  isothermal  transformation  diagram  is  shown  in 
Fig.  1  together  with  a  cooling  curve  of  M°F/sec.  The  cooling  curve  is 
constructed  to  start  from  the  (eutectoid)  temperature  if  the  trans¬ 
formation  product  is  the  eutectoid  microconstituent  (i.e.  pearlite  or 
bainite),  from  the  Ag  temperature  if  the  product  is  primary  ferrite,  or 

from  the  A  if  it  is  proeutectoid  cementite. 
cm 

The  cooling  curve  in  Fig.  1  intersects  the  curve  representing  the 
start  of  isothermal  transformation  at  a  point  X  corresponding  to  a 
temperature  T^  and  a  time  t^.  An  arbitrary  lower  temperature  Tq  (on  the 
cooling  curve)  is  chosen  which  is  reached  after  cooling  for  t  seconds. 

Two  basic  assumptions  are  proposed  by  Grange  and  Kiefer  in  regard  to 
the  transformations  occurring  between  the  points  X  and  0,  they  are: 

1.  The  extent  of  transformation  of  austenite  from  the  start  of 

cooling  to  the  temperature  T^  is  the  same  as  if  the  steel 

had  been  quenched  rapidly  from  the  austenitizing  temperature 

to  T  . 
x 

2.  On  cooling  through  the  limited  temperature  range  T  to  T  , 

X  o 

the  amount  of  transformation  is  approximately  equal  to  that 
which  would  transform  isothermally  at  the  mean  temperature 
T,.  =  4(T  +  T)in  the  time  interval  tA  =  t  -  t  „ 

Of  the  two  assumptions,  it  is  the  second  which  allows  the  graphical 
calculation  of  the  start  and  finish  curves  of  the  continuous  transformation 
diagram.  The  choice  of  a  point  0  which  corresponds  to  a  point  on  the 
start  or  finish  curve  of  the  continuous  transformation  diagram  requires 
that  the  time  tft  =  t  -  t^  be  equal  to  the  time  required  for  the  start  or 
finish  of  isothermal  transformation  at  the  mean  temperature  T.,{.  It  is 
this  which  makes  necessary  the  many  approximations  in  choosing  the  proper 
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value  of  T  ,  This,  in  turn,  renders  the  problem  suitable  for  the  high 
speed  digital  computer „  This  same  method  of  calculation  applies  to  all 
curves  represented  on  the  isothermal  transformation  diagram  with  the 
exception  of  the  martensite  start  and  finish  temperatures. 

A  basic  program  has  been  written  for  the  computation  of  the  continuous 
transformation  diagrams  from  isothermal  transformation  data.  All  that  is 
needed  as  input  data  are  the  coordinates  of  several  points  on  the  iso¬ 
thermal  curve  above  the  nose  temperature,  several  cooling  rates,  and  the 
temperature  at  which  the  cooling  curves  are  to  start.  The  isothermal 
curves  are  approximated  by  straight  line  segments  and  the  data  points 
must  be  chosen  so  that  they  are  closer  together  near  the  nose  than  at 
longer  times. 

This  program  has  been  run  with  data  for  the  SAE  1080  and  SAE  4340 
steels.  In  Fig.  2  and  3  are  presented  the  isothermal  transformation 
diagrams  and  the  resulting  continuous  transformation  diagrams  which  are 
compared  to  those  continuous  transformation  curves  determined  experimentally. 
The  agreement  between  calculated  and  experimental  diagrams  is  good  and 
bears  out  the  assumptions  used  by  Grange  and  Kiefer. 

The  basic  program  is  presented  in  the  Appendix  to  this  report  together 
with  detailed  instructions  concerning  the  selection  of  the  input  data. 

It  is  hoped  that  the  availability  of  this  program  will  lead  to  the 
greater  use  of  continuous  transformation  rather  than  isothermal  trans¬ 
formation  diagrams  in  the  heat  treatment  of  steels  and  the  consequent  im¬ 
provement  in  the  efficiency  of  this  important  operation. 
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Appendix 

Basic  Fortran  Program 


The  basic  program  for  computation  of  continuous  transformation 
diagrams  from  isothermal  data  is  presented  on  the  following  pages. 

The  program  has  been  written  in  FORTRAN  IV  and  it  is  hoped  that 
sufficient  comment  cards  were  inserted  to  make  it  understandable  to 

t  • 

anyone  familiar  with  programming . 

The  input  data  consists  of  the  following: 

(1)  Title  of  curve  to  be  calculated 

(2)  Data  for  isothermal  curve 

This  is  expressed  in  terms  of  N  data  points  ISTHER  (I,  J) 
where  I  refers  to  the  time  and  J  the  temperature  coordinate. 

(3)  The  temperature  TEUT  from  which  cooling  is  to  start 

(4)  A  set  of  cooling  curves 

This  is  expressed  in  terms  of  M  cooling  rates  of  MT  degrees 
Farenheit  per  second  (expressed  as  positive  values). 

(5)  A  suggested  step  size,  and  criteria  for  X  and  Y  points 

This  is  expressed  as  a  step  in  time  to  choose  a  value  of 

point  0  in  Fig.  1.  The  criteria  for  X  and  Y  represent  values 

of  X  and  Y  such  that  if  tA  =  t  -  t  falls  within  these  values 

**  o  x 

of  the  isothermal  curve,  the  value  is  considered  acceptable. 

(6)  A  maximum  number  of  trials  to  make  for  each  cooling  curve. 

This  has  been  25  in  the  calculations  shown  in  Figs.  2  and  3. 

The  isothermal  input  data  is  used  to  approximate  the  isothermal 
transformation  curve  into  linear  segments  and  consequently  data  points 
should  be  more  closely  spaced  near  the  nose  region  of  the  curve.  The 
cooling  curves  should  range  from  that  which  is  just  tangent  to  the  iso¬ 
thermal  curve  to  a  rate  sufficiently  slow  that  no  further  transformation 
would  be  expected,  e.g.  for  the  calculations  performed  above  6  or  7 
cooling  rates  were  employed. 


The  step  size  and  critical  values  of  X  and  Y  depend  on  the  time 

scale  of  the  curve.  If  the  curve  is  the  4340  bainite  start,  the  times 

involved  are  very  short  and  consequently  a  small  step  size,  perhaps  5 

seconds  would  be  used.  If  the  curve  were  displaced  so  the  times  were 

3  5 

of  the  order  of  10  to  10  ,  larger  steps  should  be  chosen.  The  critical 
values  of  X  and  Y  depend  on  the  position  also  and  must  be  specified 
accordingly.  The  smaller  these  values  the  more  approximations  may  be 
necessary  before  the  computed  value  is  acceptable  and  a  critical  judge¬ 
ment  must  be  made  before  these  values  are  inserted.  For  the  4340  bainite 
start,  the  critical  values  were  1  second  and  1  degree. 

As  a  further  precaution,  it  is  wise  to  place  several  blank  data 
cards  in  the  deck  to  insure  that  operation  is  terminated. 
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